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Abstract. Generalizing similar results for viscous shock and relaxation waves, we establish 
sharp pointwise Green function bounds and linearized and nonlinear stability for traveling 
wave solutions of an abstract viscous combustion model including both Majda's model and 
the full reacting compressible Navier-Stokes equations with artificial viscosity with general 
multi-species reaction and reaction-dependent equation of state, under the necessary condi- 
tions of strong spectral stability, i.e., stable point spectrum of the linearized operator about 
the wave, transversality of the profile as a connection in the traveling-wave ODE, and hyper- 
bolic stability of the associated Chapman- Jouguet (square-wave) approximation. Notably, 
our results apply to combustion waves of any type: weak or strong, detonations or defla- 
grations, reducing the study of stability to verification of a readily numerically checkable 
Evans function condition. Together with spectral results of Lyng and Zumbrun, this gives 
immediately stability of small-amplitude strong detonations in the small heat-release (i.e., 
fluid-dynamical) limit, simplifying and greatly extending previous results obtained by en- 
ergy methods by Liu-Ying and Tesei-Tan for Majda's model and the reactive Navier-Stokes 
equations, respectively. 



1. Introduction 

In this paper, we extend the viscous shock stability theory of [ZH, MaZ2, MaZ3, MaZ 1 
IHZj to traveling waves of combustion models, including the simplified combustion model of 
Majda, and an artificial viscosity version of the reacting Naver-Stokes equations. Specifically, 
we (i) derive sharp pointwise Green function bounds, yielding a sharp L x f\LP — > LP linearized 
stability criterion in terms of an Evans function condition, and (ii) assuming the Evans 
stability condition, establish nonlinear stability for waves of arbitrary type: weak or strong 
detonation, weak of strong deflagration. 

The results described in this paper, Theorems 11.21 and 11.51 below, represent in particular 
the first stability results of any kind for large-amplitude combustion waves and for weak 
detonations of Majda's model. 

This reduces the question of linear and nonlinear stability to verification of the Evans 
condition, a criterion that is readily checked numerically |Brl [ IBr2"l IBr3|. IBrZt IEDG . 

1.1. Combustion models. We show that viscous shock and combustion waves, like their 
hyperbolic counterparts, can be studied within a common framework. Indeed, viscous shocks, 
viscous detonations, and relaxation shocks may all be considered as traveling waves of the 
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special class of hyperbolic-parabolic balance laws, or reaction-diffusion-convection equa- 
tions, 



(1.1) U t + T{U) = 0, T{U) = F(U) X - (B(U)U X ) X - G(U), 
having the damping property 

(1.2) ft<r{dG)<Q, 

where (here and elsewhere) a(M) denotes spectrum of a matrix or linear operator M. For 
viscous shocks, G = 0, while for relaxation shocks, dG has constant rank, its kernel corre- 
sponding to a local equilibrium manifold. 

By contrast, combustion equations have the composite structure 

G(U) = <P(U)G(U), 

where is a scalar "ignition function" that turns the reaction on or off- specifically, it is 
zero on some subset of the state space and positive elsewhere- and G is a relaxation type 
term, dG has constant rank and a(dG) < 0: that is, an interpolation between the viscous 
and relaxation case. Thus, traveling combustion waves exhibit features of both viscous and 
relaxation shocks, in various different regimes, and our analysis must take this into account. 
Specifically, we study a subclass of (jl.lj) . comprising systems of the form, 



1.3) 



Ut + f(u, z) x = bu xx + qk(f>(u)z, 
z t = dz xx - k(f>(u)z, 



where aeK" and z G !R r , and is a "bump" -type ignition function. The physical constant 
q is the heat release parameter. Here, q > corresponds to an exothermic reaction. 

When n = 1 and r — 1, (|1.3|) is Majda's single-reaction combustion model. Then, u 
is a lumped variable combining various aspects of specific volume, particle velocity, and 
temperature, while z G [0, 1] is the mass fraction of reactant. The positive constant k 
represents the rate of the reaction. In Majda's model, the diffusion coefficients b and d are 
also assumed to be positive constants. In the following, we scale the variables so that 6 = 1. 

The vectorial version of ([1.3)1 . with u G M. n , z G lR r , and b and d positive definite matrices, 
is sufficient to encompass the artificial viscosity version of the full reactive compressible 
Navier-Stokes equations written in Lagrangian coordinates, with multi-species reaction and 
reaction-dependent equation of state, where u = (r, v, E), with r, v, and E denoting specific 
volume, velocity, and energy density, (p = <p(T), z±, . . . , z r denoting mass fractions of reactant 
species, and k matrix- valued with eigenvalues of strictly negative real part |Z1| |LyZ 2|. 

Throughout the paper, we shall carry out in parallel the analysis of the scalar and the 
(artificial viscosity) system case, exposing the main ideas in the simpler setting of Majda's 
model, then indicating by a series of brief remarks the extension to the general case. 

Physical (as opposed to artifical) diffusion terms are of form (b(u, z)u x ) x and (d(u, z)z x ) x 
with b, d matrix- valued and b semi-definite |LyZl[ LyZ2| . The diffusion coefficient b is 



commonly assumed to depend on u alone; however, like the equation of state, it properly 
depends on the make-up of the gas, hence on the mass fraction z of the reactant. See 
comments in section 11.31 about the extension of the results of this paper to such systems. 
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1.2. Statement of the results. Consider a general traveling-wave solution U(x, t) = U fa- 
st) of (jl.lj) . and the associated linearized equation about U in moving coordinates (x — st, t): 



(1.4) U t + LU = 0, LU := {F'{U) - sd x )U. 

Proposition 13.21 below, adapted from |LyZ2| , discusses the question of the existence of 
traveling waves for the scalar version of (jl.3j) ; regarding the system version, see Remark 
2. 



Definition 1.1. Let X and Y be two Banach spaces. A traveling wave U solution of (jl.lj) 
is said to be X —>■ Y linearly orbitally stable if, for any solution U of (jl.4j) with initial data 
in X, there exists a phase shift 5, such that U(-,t) approaches S(t)U'(-), in Y and as t — > oo. 

Our first main theorem is the following linearized stability criterion, precisely analogous 
to that of the viscous shock case, in terms of the Evans function D(-) associated with the 
linearized operator L about the wave, an analytic function defined on frequencies A : 9?A > 0, 
whose zeroes correspond to eigenvalues of L (see Section 1531 for further details). 

Theorem 1.2. A traveling combustion wave of (jl.3j) is L x f\L v — > LP linearly orbitally stable 
for p > 1 if and only if 

(1.5) D(-) has precisely one zero in {5iA > 0} (necessarily at A = 0). 

Theorem 11.21 is obtained as a result of detailed pointwise bounds on the Green function 
of the linearized equations about the wave, given in Proposition 17. II fresp. Remark 17.21 4 in 
the system case); see Section 1731 

Note that the spatial derivative U'(x) of the traveling- wave profile U is always a zero 
eigenfunction of L, a consequence of translation invariance of the original evolution equation. 
Thus, at a formal level, condition (jl.5j) corresponds to the statement that perturbations in 
all directions other than translation decay with time to the order of linear approximation, 
or linearized orbital stability. (At the formal level only, due to the absence of spectral gap 
between cx(L) and the origin A = 0; see |ZHj for further discussion.) 

More precisely, it was shown in |LyZl| , |LyZ2| , that, similarly as in the viscous shock case 
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(1.6) D(A)= 7 A(A)+ (|A|) 

for |A| sufficiently small, where 7 is a constant, and A is a homogeneous degree one Lopatinski 
determinant. 

The constant 7 measures the angle between the unstable subspace at —00 and the stable 
subspace at +00 for the traveling wave ODE (in this paper, equations (j3.1j) - ()3.3j) ). that is, 
the transversality of the traveling wave as a solution of the traveling wave ODE. 

The condition A (A) 7^ for 3?A > and A 7^ is equivalent to linear stability of the 
corresponding inviscid shock (square wave approximation) as a solution of the hyperbolic 
Chapman- Jouguet equations (the Chapman- Jouguet limit is the instantaneous reaction limit, 
or k = +00, of the inviscid (b, d = 0) limit of (jl.3|0 . 
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Thus, similarly as in the viscous shock or relaxation case, condition (|1.5|) is equivalent to, 



a{L) C < 0} U {0}, 
;i.7) { 7^0, 

A (A) ^ for 3?A > and A ^ 0, 

that is, strong spectral stability (first condition in ()1.7|) ). plus transversality, plus Lopatinksi 
stability of the associated square-wave ( Chapman- Jouguet) approximation. 
Note that (J1.6)) holds in the much more general multidimensional weii [zn ucwj . 



Remark 1.3. It is shown in |JLWj that, under "standard" assumptions of a reaction- 
indepenent, ideal gas equation of state, strong detonations are always Chapman- Jouguet 
stable. Together with (jl.6|) . this has the interesting consequence that transition from vis- 
cous stability to instability as parameters are varied must occur either by breakdown of 
transversality in the traveling-wave connection, or else by crossing of the imaginary axis 
of one or more nonzero complex conjugate eigenvalue pairs, i.e., a Poincare-Hopf type bi- 
furcation. This agrees with physically observed "galloping" or "pulsating" instabilities; see 
|LyZ21 ITZU ITZ2j for further discussion. 

Definition 1.4. Let X and Y be two Banach spaces. A traveling wave solution U of (jl.lj) 
is said to be X — > Y nonlinearly orbitally stable if, for any solution U of (|1.1|) with initial 
data sufficiently close in X to U, there exists a phase shift 5, such that U(-,t) approaches 
U(- — 5{tj), in Y and as t — ► oo. //, also, the phase S(t) converges to a limiting value 5(+oo), 
the profile is said to be nonlinearly phase- asymptotically orbitally stable. 

Using the information given by Theorem II .21 we further obtain our second main theorem, 
asserting that strong spectral stability implies nonlinear stability. This is a corollary of the 
more detailed, pointwise version given in Theorem 11.51 in which we let, 

L°° '■— {f G S'(M), (1 + | . |)3/2/(.) e 

In particular, L°° <^-> L H L p , for all 1 < p < +oo. 

Theorem 1.5. Under condition (|1.5j) . a traveling combustion wave U(x — st) of (I1.3J) is 

L°° L p nonlinearly phase- asymptotically orbitally stable, forp > 1. More precisely, given 
U a traveling-wave solution of (jl.3J) . given 1 < p < oo, there exist E > 0, C > 0, 8(-) G C 1 , 
and S(+oo) G R, such that the unique solution U of (II. 1JI issuing from the initial datum 
U + Uq, where, 

(l + \x\) 3 / 2 \U (x)\<E , 

satisfies the asymptotic estimates, 

\U(x,t) - U(x - 5(t))\ LP < CE (1 + tyti x -*\ 
(1.8) <CEo(l + t)-\ 

\5(t) - 5 (+oo)\ < CE (l + t)- 1 ' 2 . 



1.3. Comments. We indicate in this section how the above theorems relate to previous 
mathematical results on combustion waves. 

Strong detonations are combustion waves for which the underlying gas dynamical shock 
is of Lax type (see section l3~Tj) . 

Nonlinear stability of small-amplitude strong detonations in the small-g limit was estab- 
lished by Li, Liu and Tan |LLTj using spectral analysis together with Sattinger's method of 
weighted norms |Saj and by Liu and Ying |LYij using energy estimates, for Majda's model. 
Nonlinear stability of strong detonation in the small-g limit for the related Majda-Rosales 
model (where Zt is replaced by z x in the reaction equation), with explicit rates of conver- 
gence, was established by Li |Lilj . Nonlinear stability of small-amplitude strong detonations 
in the small-g limit for full reactive Navier-Stokes (with Heaviside-type ignition function, 
and reaction-independent equation of state) was established by Tan and Tesei in |TTj using 
detailed energy estimates. 

Roquejoffre and Vila |RVj studied spectral stability of arbitrary amplitude strong deto- 
nations in the small-fc (ZND) limit, for Majda's model (in the case d — 0). Together with 
the weighted norm argument of LLTj, this is sufficient to yield nonlinear stability with 
time-exponential rate, for exponentially decaying initial data. 

Most recently, Lyng and Zumbrun [LyZl] have shown by an elementary perturbation argu- 
ment using an abstract Evans function framework that spectral stability of strong detonations 
in the small-g (i.e., fluid-dynamical) limit amounts to spectral stability of the underlying gas- 
dynamical shock. 

For strong detonations, in the Majda model case, the spectral results of jLLTj yield (jl.5|) 
and thus full linearized and nonlinear stability. We note that, even though the weighted 
norm method suffices (as pointed out in |LLTj ) to yield a nonlinear stability result for strong 
detonations with exponentially decaying initial perturbations, our result applies to much 
more general (in particular, algebraically decaying) data and yields additional pointwise 
detail on solution structure. 

As in the shock wave case, our approach yields ultimately the reduction of stability analysis 
to a spectral problem. Thus, the following result is a consequence of Theorem 11.51 together 
with the spectral analysis of Lyng and Zumbrun LyZl|. This greatly extends and simplifies 



the strongest prior result of jTTj, illustrating the power of the method. 

Corollary 1.6. Strong detonation waves of ()1.3|) are L°° — > L p nonlinearly orbitally stable, 
p > 1, in small-q limit if and only if the limiting gas-dynamical shock is stable; in particular, 
for Majda's model, they are always stable. 

Remark 1.7. It would be interesting to extend the spectral analysis of [RV to the case 
d ^ 0, which would then imply nonlinear stability of arbitrary strength strong detonations 
for Majda's model in the small-A; (ZND) limit. More interesting still would be to extend 
this to the system case. We conjecture that the proper system analog, similar to the small-g 
result of |LyZl| , is that stability in the ZND limit is equivalent to gas-dynamical stability of 
the component Neumann shock (see discussion, |LyZl| , IU3T| IGS2j ) together with hyperbolic 
(i.e., Lopatinski) stability of the associated ZND detonation. This should be provable by a 
combination of the singular perturbation methods of |PZl IFSj and "multi-pulse" calculations 
carried out for multiple traveling-pulse solutions in models of nerve-impulse and optical 
transmission. This would recover the jRVj result of stability for the scalar, Majda model, 
for which the gas-dynamical shock, since scalar, is automatically stable (see, e.g., IS a]) and 
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the ZND detonation may readily be calculated to be stable. For systems, however, ZND 
detonations are often unstable, so that stability in the ZND limit should not be expected. 

Weak detonations are combustion waves for which the underlying gas dynamical shock is 
undercompressive (see section l3~Tj) . 

Nonlinear stability of weak detonations was obtained by Szepessy jS] for small-amplitude 
waves with intermediate k, and by Liu and Yu |LYj for arbitrary amplitude waves in the 
large- limit. Both these papers deal with the Majda-Rosales system. 

As far as we know, Theorem II .51 is the first analytical result on nonlinear stability of weak 
detonations for the Majda model (more generally, the vectorial version including reactive 
Navier-Stokes equations with artificial viscosity), and also for deflagrations of any type. 1 

Deflagrations, weak and strong, are other types of "undercompressive" combustion waves. 

It would be very interesting to determine (presumably by numerical computations) the 
existence (here assumed) and stability or instability (that is, verification of condition (|1.5jl ) 
of weak or strong deflagrations. 

Finally, note that abstract one-dimensional stability results on deflagrations are likely not 
to be so physically important, since multi-d transverse instabilities appear to play such a 
prominent role in their behavior jBj. 

As noted in [Mj . one-dimensional deflagration waves feature a pressure and a velocity 
which are locally nearly constant. Then, a "constant density" approximation |MSj shows 
that the interaction between the chemical reaction and fluid dynamics may be neglected. 
That is, roughly speaking, the complicated equations modeling reacting gas decouple into 
a part describing the fluid flow and a part describing the chemical reaction. As a result, 
deflagration waves are often modeled as systems of reaction-diffusion equations. 

The fact that detonations are usually approximated by reaction-convection equations 
(ZND), and deflagrations by reaction-diffusion equations, reflects the general belief that 
these are dominating effects in the two different contexts. Our analysis here via reaction- 
convection-diffusion puts both on the same footing, allowing treatment in a unified theo- 
retical/numerical framework, investigation/validation of these beliefs, determination of their 
realms of validity. 



1.5. Notes on the proof. An important aspect of the Lax shock analysis is that differ- 
entiated source leads to faster temporal decay |Z31 IMaZ3j . Where > c > 0, note as in 
relaxation case that source in nonequibrium mode gives faster-decaying response, at differen- 
tiated rate, and so can be treated as in relaxation case. Where = 0, undifferentiated source 
does not appear, and so can be treated as in usual shock case. What makes this technically 
feasible is that, near traveling waves, the two regions are spatially separated, corresponding 
to x < —M and x > M, M > 0, respectively. The intermediate regime Co > > is 
localized within the internal layer, corresponding undifferentiated source appears with expo- 
nentially decaying multiplier e~ e)[X \ 9 > 0. But, sources of the latter order appear already in 
the undercompressive shock case, and can be treated by the methods of |HZj with no change. 

^dn particular, we note that the weighted norm technique of jLLT] is inherently restricted to strong 
detonations for the scalar, Majda model; see the discussion of TZ1 . 
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So, our analysis is by interpolation between the viscous undercompressive shock analysis of 
|HZ| IMaZ3j and the relaxation shock analysis of jMaZlj . The new aspects of the argument 
not present in the undercompressive shock case are isolated to bound (|8.8|) (ii). Remark 
18.21 and the new auxiliary Lemma 18.61 In particular, no new convolution estimates were 
necessary, only a series of observations having to do with the fact that undifferentiated sources 
appear in a direction for which the Green function decays more rapidly, at differentiated rate. 

Remarks 1.8. 1. As a relaxation system, (J1.3|) is somewhat degenerate, violating the 
usual assumption of genuine coupling between equilibrium and relaxation variables u and z 
associated with time-asymptotic smoothing of solutions (see, e.g., jMaZl|IZ2] and references 
therein). Indeed, asymptotic decoupling of the ^-equation plays an important role in the 
analysis; see Remarks 14 . 1 1 and 14 . 41 2 . Diffusion terms b, d not present in standard relaxation 
systems enforce smoothing directly. 

2. The case d = that is often considered for Majda's model requires slightly different 
handling. Absence of ^-diffusion leads to "hyperbolic" delta-function components in reac- 
tive modes reminiscent of those encountered in |MaZl| IMaZ3j in the case of relaxation or 
degenerate viscosity, but with the difference that incoming modes on side x > are not 
time-exponentially damped. This can be accomodated in the analysis by the introduction 
of an exponentially weighted norm in the spirit of |Saj in the ^-component only, for x > 0, 
using the property of exponential decoupling as x — > of reactive and fluid modes. This is 
a mathematical issue only; for physical models, d is strictly parabolic: 3la(d) > 0. 

1.6. Extension to the Navier-Stokes equations with physical viscosity. The full 

reactive Navier-Stokes equations with real, or physical viscosity may be treated by essentially 
the same techniques, using the more complicated arguments (and more detailed Green fn. 
bounds) developed in |MaZ31 IMaZ41 IZ2l El IHRZ] for the treatment of viscous shocks with 
real viscosity. However, these arguments so far are limited to the strong detonation case. 
(Likewise, for technical reasons, the viscous shock theory is so far limited for physical viscosity 
to the Lax and overcompressive case.) We leave this to a future work. 

1.7. Plan of the paper. In Section El we describe Majda's combustion model and its 
vectorial generalization, in Section 01 the various types of traveling wave connections that 
may occur, and in Section 0] the linearized eigenvalue equations about these traveling waves. 
In Sections and |H1 we construct the Evans function and resolvent kernel of the Linearized 
operator about the wave following the abstract framework of ZH, MaZ3 , specializing in the 
low-frequency regime to the special structure of (jl.3|) using the limiting constant-coefficient 
calculations of Section 0J In Section [?1 we convert the resulting pointwise resolvent kernel 
bounds to pointwise Green function bounds by stationary phase type estimates on the Inverse 
Laplace transform formula, in the process establishing Theorem II .21 equating linearized and 
spectral (Evans) stability. Finally, in Section |H1 we carry out a nonlinear stability analysis, 
establishing Theorem 11.51 

2. Preliminaries 

2.1. Majda's model. We begin with the scalar version of system (jl.3|) . We assume as in 
|LyZ2l that /, <p G C 2 , 

(2.1) f u (u,z)>0, f uu (u,z)>0, 
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and that is a "bump" -type ignition function that is identically zero for u < Ui or u > u l 
and strictly positive for m < u < u 1 . 

It is sometimes useful to rewrite (|1.3|) in the conservative form 



(2.2) (u + qz) t + f(u, z) x = u xx + qdz x 

(2.3) z t = dz xx - k<p(u)z. 

Remarks 2.1. 1. Note that the flux /, modeling equation of state, depends on z, model- 
ing the chemical constituation of the gas, with the linear dependence loosely following the 
averaged equation of state derived in jCHTj for the full Euler equations. This is important 
for realistic modeling of the full equations of reacting flow; see [CUT, LyZ2 for further dis- 
cussion. For the Majda model, new qualitative phenomena emerge for f z ^ at the levels 
of both existence and behavior of detonation profiles |LyZ2| . 

2. Following [MJ, <fi is usually taken to be a "step" -type function, vanishing for u <ui and 
positive for u > ui. As discussed in |LyZ2| , our alternative choice of a bump-type function 
is motivated by the physical parametrization of temperature with respect to velocity u in 
the traveling- wave phase portrait of the ZND model. This choice admits all the phemomena 
of the step-ignition case, restricting to u < u l . In addition, it allows for existence of weak 
deflagration profiles (defined in Section EJ), as the step- type ignition function in general does 
not; see |LyZ2| or Remark 13.11 

2.2. Reacting Navier— Stokes equations. The single-species reacting Navier-Stokes equa- 
tions with artificial viscosity, written in Lagrangian coordinates, take the form 

Tt~V x = 0, 

v t +Px = hv xx 

(2 4) 

E t + (pv) x = b 2 E xx + q 3 k(f)(T)z, 
z t = dz xx - k(f)(T)z, 

where r, v, E = e+qz+v 2 /2, z denote specific volume (p _1 , where p is density), velocity, total 
energy density, and mass fraction of reactant, T = T(r, e) temperature, and p = p(r, e, z) 
pressure, k,q 3 ,bj,d > constant, or 












: 









6i 








V° 








in (II. 3|) . The ignition function is assumed to vanish identically for T below some critical 
ignition temperature Tj, and to be strictly positive for T above Tj. 

A common choice of equation of state is the reaction-independent gamma-law 

(2.5) p = Te/r, T = c _1 e, 

where c is the specific heat constant and T = 7 — 1 is the Gruneisen constant. In the 
thermodynamical rarified gas approximation, 7 > 1 is the average over constituent particles 
of 7 = (n + 2) /to, where to is the number of internal degrees of freedom of an individual 
particle |Baj . 
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A more accurate assumption, following |CHTj . is to view the g composite of un- 

burned and burned phases with different equations of state Tjfa, ej), Pj(rj, ej), j = 1 corre- 
sponding to the unburned and j = 2 to the burned state, with 

7i = r/z, t 2 = r/(l - z) (i.e., p x = zp, p 2 = (1 - z)p), 

T\ =T 2 = T, e = ei + e 2 , p = Pi + P2- 

If both phases obey (different) gamma-law equations of state, this leads |CHTj to a gamma- 
law-type equation of state ()2.5|) with reaction-dependent coefficients 

(2.7) c(z) = «* + (1 - z)c 2 , r W = !££i±lizLil££i . 

ZCi + (1 - Z)C 2 

This is the "typical" equation of state we have in mind. 

3. Traveling Waves 
We consider traveling- wave solutions, i.e., solutions of the form 

u(x, t) = u(x — st), z(x, t) = z(x — st), s > 0, 

of ()1.3|) that connect an unburned state (u + ,z + ) = (u + , 1) to a completely burned state 
(m_,2_) = (w_,0). These are combustion waves that move from left to right leaving com- 
pletely burned gas in their wake. 

Thus, dropping bars for notational convenience, we find that the traveling-wave Ansatz 
leads, after an integration, from (jl.3j) to the system of ordinary differential equations: 

(3.1) u = f(u, z) — /(m_, Z-) — qdy — sqz — s(u — w_), 

(3.2) z! = y, 

(3.3) y' = d' 1 ^ — sy + k(p(u)z), 

where we have used y := z' to write the system in first order and ' denotes differentiation 
with respect to the variable £ := x — st. We assume that the end states are such that 

(3.4) m_ G [u h u% u + (£ [ui,u% 
so that 

(3.5) <f>(uJ) > 0, <f>(u + ) = 0, <f)'(u + ) = 0. 

Equation ()3.4|) has the physical interpretation that the unburned end state is below ignition 
temperature so that the there is no chemical reaction on the unburnt side. 

3.1. Rankine— Hugoniot conditions. A necessary condition for the existence of a connec- 
tion is that the end states at ±oo be rest points of the traveling-wave equation. This leads 
to the Rankine- Hugoniot condition 

(RH) f(u + , z + ) - /(u_, z-) = sq + s(u + - u_), 

together with the requirements that 

(3.6) y± = 4 = 
and (justifying assumptions ()3.5j) ) 

(3.7) <P{u ± )z ± = 0. 
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Restricting now to Majda's model, write a± := f u (u±,z±) and (3± := f z (u±,z±). Then, 
the traveling-wave profile is said to be a strong detonation if 

(3.8) a_ > s > d+. 
It is said to be a weak detonation if 

(3.9) s > a_, a+. 
Similarly, it is said to be a weak deflagration if 

(3.10) &+ > s. 
It is said to be a strong deflagration if 

(3.11) a + > s > a_. 

Degenerate profiles for which the inequalities are nonstrict are called Chapman- J ouguet 
detonations or deflagrations and lie on the boundary between weak and strong branches. 

Remark 3.1. For / independent of z, we find by (J2.1j) that detonations correspond to case 
U- > Ui > u + , deflagrations to case w_ < u % < u + . In particular, deflagrations cannot occur 
for a step-type ignition function, for which u % = +oo. 

A routine modification of Propositions 2.1 and 2. 2, |LyZ2| , accounting for ^-dependence 
of /, yields the following description of solutions of (jRH|) . 

Proposition 3.2. For fixed u + , suppose that f(u + ,l) < f(u + + q,0). Then, there exist 
s*(u + ) < s*(w + ) such that (i) for s > s* there exist two states w_ > u + for which (jRH|) 
(but not necessarily (|3.5p ) is satisfied (weak and strong detonation), for s = s* there exists 
one (Chapman- Jouguet detonation), and for s < s*, there exist none no solutions w_ > u+. 
(ii) For s < s* , there exist two states w_ < u + for which f)RH|) is satisfied (weak and strong 
deflagration), for s = s* , there exists one (Chapman- Jouguet deflagration), and for s > s* , 
there exist none. If f(u + , 1) > f(u + + q,0), on the other hand, then for each s there exists 
at most one solution w_ > u + and one solution w_ < u + , (strong detonation and strong 
deflagration, respectively). 2 

Remark 3.3. The case /(«+, 1) < f(u+ + q, 0) is essentially identical to that of the reaction- 
independent case discussed in |LyZ2 , for which f(u + ) < f(u + + q) by ()2.1j) . More generally, 
f z < is sufficient for f(u + , 1) < f(u + + q, 0). For further discussion of the (jRH|) problem, 
see |LyZ2| : in particular the Chap man- Jouguet diagrams of Figure 1. 

3.2. The connection problem. Linearizing (|3.1|) - (|3.3JI around the state (w_, z~, y~), we 
find the constant-coefficient system of ordinary differential equations 

(u\ /&- — s — sq —qd \ ( u\ 

z) = 1 )iz . 

y) \ kd- x cj>{uJ) -sd- 1 ) \yj 

For strong detonations, the coefficient matrix in ()3.12j) is easily seen to have two positive 
eigenvalues and one negative eigenvalue. Thus, there is a two-dimensional unstable manifold 



Typically, both; in particular, if / grows superlinearly in |it|, then both solutions exist for each s. 
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at (w_,0,?/_). Similarly, we note that the linearization of ()3.1|) - ([3.3|) about (u + , z + ,y + ) is 

(u\ I a+ — s b + — sq —qd \ (u 
z = 1 \\ z 

y) \ o -sd- 1 ) \y 

By the block-triangular structure, it is easy to see that there are two negative eigenvalues 
and one zero eigenvalue. It is easy to see that the center manifold is a line of equilibria, so 
plays no orbit may approach the rest point (u + , z + ,y + ) along the center manifold. Thus, for 
connections, the important structure is the two-dimensional stable manifold at (u + ,l,y + ). 
Counting dimensions, we see that a connection corresponds to the intersection of two two- 
dimensional manifolds in R 3 . In particular, it generically persists as a unique, transverse 
intersection, under variations in parameters such as «±, s consistent with (jRH|) . See |LyZ2 
for a discussion of this situation in the case d = 0. 

Similarly, for weak detonations, there is a one-dimensional unstable manifold at (u_, 0, yJ) 
and a two-dimensional stable manifold at (u + ,0,y + ). Thus, connections are typically codi- 
mension one in the set of (jRHJ) compatible parameters, in contrast to the strong detonation 
case. See |LyZ2| in the case d = 0. This situation is analogous to that of a Lax-type shock in 
the nonreactive case. For weak deflagrations, there is a a two-dimensional unstable manifold 
at (m_,0,?/_) and a one-dimensional stable manifold at (u + ,0,y + ). Thus, connections are 
again generically codimension one in the set of (jRH|) compatible parameters. See [LyZ2 
in the case d = 0. For strong deflagrations, there is a a one-dimensional unstable mani- 
fold at (u-,0,y-) and a one- dimensional stable manifold at (u + ,0,y + ), and connections are 
codimension two. 

In every case, we have by the discussion above: 

Lemma 3.4. Traveling-wave profiles (u, z) corresponding to weak or strong detonations or 
deflagrations satisfy 

(3.14) \(d/dx) k ((u,z) - (u,z)±)\ < Ce~ 9lx \ x^O, < k < 3. 

Proof. Standard ODE estimates for stable and unstable manifolds. □ 

Remark 3.5. Weak detonations and deflagrations are analogous to undercompressive 
shocks in the nonreactive case, with strong deflagrations undercompressive of degree two; see 
|Z2j for a discussion of shock classification. In the case d — 0, it can be demonstrated that 
weak detonation connections do occur in some cases, but deflagration connections (weak or 
strong) of the type we have described do not |LyZ2| . 3 It is an interesting question whether 
or not they occur for d ^ 0, or, more generally, for the full, reactive Navier-Stokes equations 

|LyZTl|L^2l . 

3.3. The system case. Under the further assumption of (asymptotic) dissipativity, 

(3.15) 9fe7(di£ - eb)± < 4t := (df/du)(U±), 

i + Kr 

9 > 0, for all £ G R (standard for systems |Z2j). it is readily verified using the above- 
mentioned block-triangular decomposition of limiting systems into fluid and reactive blocks 

3 More precisely, weak deflagration profiles exist only in the degenerate case u + = u' ', with u(x) converging 
to ?i_ as x — > +oo at subalgebraic rate; strong deflagration profiles do not exist in any case |LyZ2| . 
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together with fluid dynamical results of Majda and Pego |MPj , that the main results of this 
section carry over to the system case, substituting for definitions (|3.8|) - (|3.11|) the system 
versions (for right-going waves, s > 0) 

(3.16) a~ > s > a£, s > af , j ^ n, (strong detonation) 

(3.17) s > a~, &n, s > Sty, j ^ n, (weak detonation) 

(3.18) a~,a^ > s, s>af,j^n, (weak deflagration) 
and 



(3.19) > s > a n . s > , j ^ n, (strong deflagration) 

where af < ■ ■ ■ < denote the eigenvalues of 6c± := f u (u±, z±). 

In particular, connections if they exist satisfy ()RH|) . and if they are transverse are gener- 
ically of codimension zero, one, one, two, respectively, in the set of (|RH|) -compatible pa- 
rameters. Further, the connecting profile satisfies (|3.14|) . converging exponentially to its 
endstates t/± as x — > ±oo. See |Zlj . Appendix A for further details. Condition (|3.15|) holds 
trivially for identity viscosity b — I, and holds also for the physical (semidefinite) viscosity 
of the reacting Navier-Stokes equations [MaZlJ IZ2j . the two main cases we have in mind. 
For simplicity of notation, we assume also that a(d) is semisimple, so that d has a full set of 
eigenvectors. 

Likewise, there is a simple analogy to the Chapman- Jouget analysis of Section 13.11 and, 
for the typical mixed gamma-law equation of state ()2.7|) of Proposition 13.21 For, rearranging 
flRH|) in the case of the reacting Navier-Stokes equations (J2.4)) . we obtain |LyZl| from the 
first equation that (v+ — vJ) = — s(r + — r_), from the second that (r,p)± lie on the Rayleigh 
line 

(3.20) p+ -p_ = -s 2 (T + -T_), 
and from the third the shifted Hugoniot curve 

(3.21) (e+ - e_) + q = (-l/2)(p+ + p^)(r + - r_). 

Thus, fixing (t,v,E) + , viewing (|3.21|) as determining a "burned" pressure law 

(3.22) p _=p_( T _), 

and assuming that e + can be recovered from r_, p_ through inversion of p_ = p(r_, e_, 0), we 
find that the allowable states (r, v, E)_ are determined as the intersection in the r — p plane 
of line (|3.20j) with curve (J3.22)) . similarly as in the scalar (Majda's model) case. Moreover, 
there are two distinct solution structures, according as 

(3.23) p + < P_(r + ) 

(standard: for P_ convex, pairs of weak/strong detonations, deflagrations as in the scalar 
case) or the reverse (nonstandard: for P_ convex, single strong detonation, deflagration). 
For the typical equation of state (j2.5|) - (|2.7j) . it is readily calculated that 

{r + /T l -{l/2){r-T + ))p + + q 



(3.24) P_(r) 



r/r 2 -(l/2)(r 
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whence ()3.23|) reduces to 



(3.25) p+(l - r 2 /ri) < qT 2 /r + . 

From (|3.25p . we see that T\ < T 2 (in particular, including the reaction-independent case) 
implies a standard (jRH|) solution structure. Roughly speaking, this corresponds to a reaction 
in which complicated compounds break up into simpler components, so that n decreases 
and r = 2/n increases, the reverse situation to a reaction in which simple components 
combine into more complicated molecules. We conjecture, by analogy with the scalar case, 
that for (|2.7|) . condition (|3.25|) equivalent to T z < implies further a standard connection 
structure, at least in the ZND limit (for which there is a close connection to Majda's model 
LyZT, LyZ2t IUST] ). However, even for I\ > T 2 , the standard (|RH|) solution structure is 



recovered for q sufficiently large. 

Remarks 3.6. 1. Dissipativity, (j3.15J) . implies in particular hyperbolicity of the first-order 
convection terms, i.e., that a has real, semisimple eigenvalues. We use this freely below. 

2. Existence of detonation connections for the full, reacting Navier-Stokes equations 
has been studied by Gasser and Szmolyan [GSl, GS2J using geometric singular perturbation 
techniques in the ZND (b, d — > 0) limit, and by Gardner [G] using Conley index methods. 

4. The eigenvalue equation 

Suppose (u(x — st), z(x — st)) is a traveling- wave profile of ()1.3|) as described above. We 
now begin to investigate the stability of such an object. The linearized equations about 
(u, z) in moving coordinates x = x — st, are, dropping tildes, 

(4.1) m - q(k(f>'(u)uz + k<p(u)z) + (au) x + ((3z) x = u xx , 

(4.2) z t — sz x = —kcf)'(u)uz — k<f)(u)z + dz xx , 

where a := f u (u,z) — s, (3 := f z (u,z), and u, z now denote perturbations. The eigenvalue 
equations corresponding to this linear system are thus 

(4.3) u" = (Mi - q(k(f)'{u)uz + k<p(u)z) + (au)' + (f3z)') , 

(4.4) z" = d~ l (Xz - sz + k(f)'{u)uz + k(j)(u)z) , 

Alternatively, upon substituting dz" — Xz + sz' = k<fi'(u)uz + k<fi(u)z from (J4.4j) into (J4.3|) . 
we can rewrite (|4.3|) as 

(4.5) u" = (X(u + qz) - sqz' - qdz" + (au)' + (f3z)') . 

Compare this with the remark at the end of Section We write f)4.3j) — ()4.4j) as a first-order 
system. To do so, we define W := (w, z, u', z') tT , so that (|4.3|) - (|4.4|) becomes 

(4.6) W' = A(x,X)W 
where the coefficient matrix is 

/ 

(4.7) A(x,A)= X + a' -qk(j)'(u)z (? - qk<f>(u) 

y dr x k$'(u)z d~ l X + d~ l k(j>(u) 

13 



1 





\ 





1 




a 


P 







—sd 


V 



System ()4.6j) has a limiting constant-coefficient structure, i.e., the coefficient matrix has 
limits as x —* ±00. That is, 

A(x, A) — > A±(A) as x —>■ ±00, 

and the limiting matrices are given by 

/0 1 \ 

1 



(4.8) 



A±(A) 



A —kqcj){u±) a± f}± 
\0 d~ 1 {\ + k<p{u±)) sd~ l ) 

Remark 4.1. Here, we coordinatize as W = (u, z,u' , z') tr following the general scheme of 
ZH , IMaZ3j as this is what we shall need below to establish the pointwise bounds. However, 
it is sometimes helpful to see the fluid/reaction structure in the system. To see this at the 
level of the eigenvalue ODEs, we write 

(4.9) W := (u,u',z,z') tr , 

separating the fluid (u, u') and reaction (z, z') quantities. In this labeling scheme, the eigen- 
value ODE (|4.6j) becomes 

(4.10) W' = k(x,X)W 

where the matrix A can easily be obtained from ()4.7|) by appropriately swapping entries, and 



(4.11) 



A±(A) 



/o 1 \ 

A a± —kq(f>(u±) j3± 

1 

\0 d~\\ + k<p(u±)) -sd- 1 ) 



In particular, the upper block-triangular structure of A± will be useful in the constant- 
coefficient analysis below. 

4.1. Constant— coefficient analysis. We now examine the constant-coefficient limiting 
systems W = A±(X)W (or equivalently W' = A±(X)W) in some detail. From the up- 
per block-triangular structure in (|4.11j) . it is quite straightforward to compute eigenvalues; 
they are simply the eigenvalues of the diagonal blocks. From the upper left-hand "fluid" 
block, we obtain the fluid eigenvalues 

, a ± ± Ja\ + 4A 
»f = o > 



-sd" 1 ± Vs 2 d~ 2 + Adr^X 



(4.12) 

while the lower right-hand "reaction" block contributes eigenvalues of form 

(4.13) £ 
and 

(4.14) fi~ = — : — — — from 



from A 1 



-sd- 1 ± yjs 2 d- 2 + 4(d-U + d- l k(f){u_)) 



The corresponding eigenvectors also have structure inherited from the block-triangular na- 
ture of the limiting matrices A±. In particular, as long as the fluid and reaction eigenvalues 
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remain distinct — as our calculations below show they are for small A, the corresponding 
eigenvectors take the form 

/ 

and vt = 



(4.15) 

where 
(4.16) 
and 
(4.17) 



/1\ 





V°7 

M + (A) := 



± 



v 



f4 



J 



o 1 

A a+ 



-til 



-l 



M_(A) := 



1 

A a_ 





/3 4 







We also record here Taylor series expansions of those eigenvalues of the limiting systems 
which become zero at A = 0, the so-called slow modes. These are 



(4.18) 

and 

(4.19) 



^(A) = 1a - ^a 2 + 



M ±(A)=±— A T ^A 2 + 



with associated (right) eigenvectors, now written in (u, z, u', z') coordinates, 



(4.20) 

and 

(4.21) 

where limiting fluid modes 
(4.22) 



>4W 



vf(X) 







+ 



+ ■ 



Rf 



(4.23) 



have vanishing ^-component. 

We shall also have need of the adjoint eigenvalue equation 

u" = X*u — qk(j)'(u)zu + k(j)'{u)zz — au , 

z" = dT 1 (A* - qk(j)(u)u + k(f)(u)z - f3u' + sz') 

associated with ()4.3|) - ()4.4|) . where A* denotes complex conjugate. Writing as a first-order 
system W' = A(x, \*)W, W = (u, z, u', z') tv , and studying the eigenvalues fif and eigenvec- 
tors ti- of the limiting, coefficient-matrix A_|-(A*), we have by duality that Jx 
a brief calculation yields Taylor expansion 



(pf)*, while 



(4.24) 



vf(\) 



L 



+ 



+ 
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at A = 0, where limiting left fluid modes = c(l, q) tv on the minus infinity (reactive) side 
satisfy 

(4.25) L- 1 (-q, If. 

(See Remark 14.41 1 below for further discussion, and extension to the system case.) 

Definition 4.2. The domain of consistent splitting for an ODE W = A(x, X)W with 
asymptotically constant coefficients is the set of A G C such that 

(i) the limiting matrices A-t(A) are hyperbolic, i.e., they have no center subspace, and 

(ii) the dimensions of the stable (unstable) subspaces S + (X) and ^"(A) (£/ + (A) and 
U~(X)) are the same. 

Lemma 4.3. The set {A G C : 3?A > 0} is a subset of the domain of consistent splitting. 
Proof. Immediate by inspection of formulas (j4.12)) - (j4.14|) . □ 
In fact, more can be said. The limiting matrices A± fail to be hyperbolic if 

det(A±(A) -if) = 

for some f G R. But, again using the upper block-triangular structure of A±, it is clear 
that this determinant vanishes if and only if the determinant of one of the diagonal blocks 
vanishes. This leads to the following four dispersion curves in the complex A-plane, 

\t(0 = d(-e + isd- 1 ) 1 

A;(f) = d(-i 2 + isd- 1 ) - (j){u-). 

These curves, parabolae opening into the left complex half plane, can be used to describe the 
boundary of the domain of consistent splitting. In particular, we note that Lemma \4. 31 gives 
that the open right half complex plane is contained in the domain of consistent splitting, 
and, varying A from right to left, we find that A cannot leave the domain unless we cross one 
of these four dispersion curves. Thus, the component of the domain of consistent splitting 
which contains +oo contains a set of form 

fi„ := {A : KA > maxj-^pAI^ISAI 2 }}, > 0. 

See ()6.1(J|) and Lemma f6.2l below. 

4.2. The system case. The above calculations extend in straightforward fashion to the 
system case, substituting block matrix computations for the scalar computations above. In 
particular, the block-triangular structure of (j4.11|) is maintained, reducing the computation 
of constant-coefficient modes to a computation on the upper lefthand diagonal fluid block 
that is exactly the viscous shock computation done in |ZH|. IMaZ3] . and a computation on 
the lower righthand diagonal reaction block that on the plus infinity side is a particularly 
simple (scalar convection) version of the same viscous shock computation and on the minus 
infinity side consists of fast modes that need not be resolved. 

Specifically, the Taylor expansions of slow modes at A = become 

1 2dt , 

(4-26) /,+ A = -A - — |-A + • ■ ■ , 

s s 6 
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1, . . . , m, z G 



(4.27) ^(A) :-A:-^. 

(4-28) (5*) +"' = fX 



and 

(4.29) t - / ± (A)=(^)+..., R%=( r <f 

(4.30) «?(A)=( L « )+-.., ^=c(|), L ]j= c(J/^ 

where and r r) j are left and right eigenvectors of d (now matrix- valued) and a^, L -, r^- 
are the eigenvalues and left and right eigenvectors of df /du(U±), bf = (lf,jbrfj)± (strictly 
positive, by dissipativity assumption (|H.15jl ). and d^ = (l r idr r> i) + . Note that we again have 
vanishing of the ^-component of asymptotic fluid modes Rfj, as well as the key orthogonality 
relation 

(4.31) L~ 1 >~ q 



Remarks 4.4. 1. The structural relations ()4.22jl . (|4.25jl and their system analogs (|4.29jl - 
1)4.31)1 for the asymptotic modes, play an important role in our analysis; see Remarks 16.91 1. 
18.71 and the proofs of Proposition 17.11 Lemma 18.6) and Proposition 18.11 Indeed, this is 
essentially the only structure that we use, other than the existence of Taylor expansions 
of slow modes at A = 0. These may be verified easily by substituting into the limiting, 
constant-coefficient eigenvalue systems 

u" = (Xu — qk(p±z + a±u + /3±z'), 

(A 32) 

z" = d- l (\z-sz' + k<P ± z) 

and their dual, adjoint systems the Ansatze U = e^ x R, U = e fJ, * x L, respectively, to obtain 
characteristic equations 

fi 2 — a±ji — XI n qk(f>± — (3±fi 

dfi 2 — XI r + sfi — k<fi±] 

and 

/ fi 2 + a±fi — XI n qkcj)± + f3±n 

dfi 2 — XI r + s/j, — k</>±) 

respectively, which, setting A = \i = to obtain slow mode behavior at A = 0, reduce to 

from which ()4.29j) - ()4.31)) are evident. 

2. The containment Q v C A likewise carries over to the system case, by dissipativity 
assumption ()3. 15|) and Lemma f6.2l below. 
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5. The Evans function 

We now construct the Evans function following the abstract framework of |MaZ3[IZ2"] . For 
historical origins of the Evans function, see |AGJ[ IPW| IGZj and references therein. 

5.1. The Conjugation Lemma. We first recall a central result connecting variable- and 
constant-coefficient ODE. Consider a general family of first-order ODE 

(5.1) W - A(x,X)W = 

of the form ()4.6|) . indexed by a spectral parameter A £ Q C C, where W £ C , x £ K. and 
'"" denotes d/dx, assuming (cf. Lemma E3J) 

(hO) Coefficient A(-, A), considered as a function from Q into C°(x) is analytic in A. More- 
over, A(-,A) approaches exponentially to limits A± as x — > ±oo, with uniform exponential 
decay estimates 

(5.2) \(d/dx) k (A - A±)| < C x e- e ^ tc \ for x ^ 0, < k < K, 

Cj, 9 > 0, on compact subsets of fl. 

The following asymptotic ODE result generalizes the Gap Lemma of |GZ| : for a proof, 
see, e.g.. jMZTl IMaZ31 IZ2] . 

Proposition 5.1 (The Conjugation Lemma [MZlJ). Given (hO), there exist locally to any 
given Ao £ ^ invertible linear transformations P + (x,X) = I + + (x, A) and P-(x, A) = 
/ + ©_(x, A) defined on x > and x < 0, respectively, &± analytic in A as functions from Q 
to C°[0, ±oo), such that: 

(i) For any fixed < 9 < 9 and < k < K + I, j > ; 

(5.3) \(d/d\y(d/dx) k e±\ < C^d^e-^^ 2 for x ^ 0. 

(ii) The change of coordinates W =: P±Z reduces ()5.1|) to 

(5.4) Z' - A±Z = for x ^ 0. 

5.2. Normal modes. Using Proposition 15.11 we next construct normal modes for ()5.1|) . 
Recall the domain of consistent splitting defined in Definition 14 .2[ Section 14.11 

Lemma 5.2. On any simply connected subset of the domain of consistent splitting A, there 
exist analytic bases {Vi, . . . , Vk}^ 1 and {V^+i, . . . , Vn}^ 1 for the subspaces S± and U± defined 
in Definitionl 



Proof By spectral separation of U±, S±, the associated (group) eigenprojections are analytic. 
The existence of analytic bases then follows by a standard result of Kato; see |Katj . pp. 99- 
102. ' □ 

By Lemma f5. 11 on the domain of consistent splitting, the subspaces 

(5.5) S + = span{W+, . . . , W+} := span{P + y+, . . . , P + V+} 
and 

(5.6) U~ := span{W fc - +1 , . . . , Wj} := span{P_V k ~ +1 , ...,P-Vj} 

uniquely determine the stable manifold as x — > +oo and the unstable manifold as x — * — oo of 
(|5.1|) . defined as the manifolds of solutions decaying as x —> ±oo, respectively, independent 
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of the choice of P±. More generally, Wj := (PVj)±, j = 1, . . . , N are called normal modes 
for (|Ojl . 

In the context of Majda's model (more generally, the system analog HERE), Va are 



comprised of the vectors described in Section 14.11 (resp. I4.2J1 . of which the slow modes, 
defined as those approaching the center subspace of A± as A — > 0, are t/jr, vf (resp. vfp 
v^Tj). As /asi growing and decaying modes, defined as those approaching the stable and 
unstable subspace of A±, hence spectrally separated both from each other and from slow 
modes, may always be chosen analytically in a neighborhood of A = by the same argument 
used in Lemma 15.21 we obtain by our asymptotic description of slow modes the following 
important extension. 

Lemma 5.3. For Majda's model (more generally, the vectorial version including reactive 
Navier-Stokes equations with artificial viscosity), normal modes extend to A U 5(0, r) for 
r > sufficiently small, with low-frequency (slow mode) asymptotics as described in Section 
gH (resp. [PU . 

Proof. The explicit Taylor expansions of Section |4~T1 (resp. 14. 2 j) yield analytic extensions on 
.8(0, r) of slow modes V^, spanning invariant subspaces of A±. As fast modes always have 
such analytic extensions, we may combine them to obtain analytic bases Vj for invariant 
subspaces S , U of A± extending those of Lemma l5~2l These in turn determine A±-invariant 
projections onto those subspaces, which must therefore be the unique analytic extension of 
the corresponding projections on An -6(0, r), and thus an analytic extension onto AUP(0, r). 
The result then follows again by the result of Kato as in the proof of Lemma 15.21 □ 

Remark 5.4. The extension of normal modes through the essential spectrum boundary to 
a neighborhood of A = is crucial for all that follows; see |GZ| IZHj for further discussion. 

5.3. Construction of the Evans function. 

Definition 5.5. On any simply connected subset of the domain of consistent splitting, let 
V-f, . . . , V k and V k+1 , . . . , be analytic bases for S + and [/_, as described in Lemma I5~2l 
Then, the Evans function for (j5.1j) associated with this choice of limiting bases is defined as 



D(X) := det (W+, . . . , W+, W k+V . . . , W N 

(5.7) ; 

= det (P+V+, . . . , P + V k + , P-V k - +1 , . . . , P_ V h 



\x=0,X 

where P± are the transformations described in Lemma f5. 11 



|x=0,A 



Remark 5.6. Note that D is independent of the choice of P±; for, by uniqueness of sta- 
ble/unstable manifolds, the exterior products (minors) P+V± A ■ • • A P+V k + and P_V k+l A 
■ • ■ A P-Vft are uniquely determined by their behavior as x — » +oo, — oo, respectively. 

Proposition 5.7 (MaZ3, Z). Both the Evans function and the stable/unstable subspaces 
S + and U~ are analytic on the entire simply connected subset of the domain of consistent 
splitting on which they are defined. Moreover, for A within this region, equation (|5.1|) admits 
a nontrivial solution W G L 2 (x) if and only if D(X) = 0. 

Remark 5.8. In the case that (|5.1|) describes an eigenvalue equation associated with an 
ordinary differential operator L, A G C 1 , Proposition 15 . 71 implies that eigenvalues of L agree 
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in location with zeroes of D. In |GJ1| IUJ2j . Gardner and Jones have shown that they agree 
also in multiplicity; see also Lemma 6.1, |ZHj . or Proposition 6.15 of |MaZ3j . 

By Lemma 15.31 we have immediately that D extends to a neighborhood of the origin. 

Lemma 5.9. For Majda's model (more generally, the vectorial version including reactive 
Navier-Stokes equations with artificial viscosity), D extends analytically to A U 5(0, r) for 
r > sufficiently small. 

6. The resolvent kernel 

Next, we estimate the resolvent kernel G\(x,y) := [L — X)~ x b y {x) associated with the 
linearized operator L about the wave, again following the abstract framework developed in 
ZH ( MaZ3j, Rewriting (jl.3|) in vectorial form 

(6.1) U t + F(U) x + G(U) = BU XX} 

(6 'u - («, W) - ( /( \f ~) , G(U) = , fl = (J 2) , 

in coordinates moving with a given traveling- wave profile U (x) = (u, z) (x) (stationary, in 
the moving coordinate frame) and linearizing about U, we obtain the linearized equations 

(6.3) U t = LU := BU XX - (AU) X + CU, 
where 

A == mm = (« _° s ) , B = (J 2) , c ^ mu) = !^,) • 

The results of |ZH| IMaZ3] for general strictly parabolic systems of the form (|6.3p . (/ 
state that the resolvent kernel is a meromorphic function on the domain of consistent splitting 
defined in Section E~T1 where it is determined by 

(6.4) (L-X)Gx = 5 y (x) 

and the property of decay as x, y — >• ±oo. Moreover, they give an explicit description of G\ 
in terms of the normal modes of the eigenvalue equation constructed in Section 15.21 from 
which we may extract sharp pointwise bounds. We cite here the relevant theory, referring 
the reader to |ZH| IMaZ3j for proof. 

6.1. Duality relation. Consider solutions U of the eigenvalue equation (L — X)U = and 
solutions U of its adjoint (L* — X)U, where 

(6.5) L*U := B tr U xx + A tT U x + C tT U 

denotes the L 2 adjoint of L and A the complex conjugate of A. Introducing the phase- 
variables W := (U, U') and W := (U, U'), write the eigenvalue equation of L and its adjoint 
as first-order ODE of form (|5.1J1 in W and W. Then, we have the following key relation. 

20 



Lemma 6.1 ( |ZHllMaZ3| ). W = (U,U') satisfies (jSHJ) if and only if 
(6.6) W*SW = constant 

for all W = (U, U') satisfying the adjoint eigenvalue equation, and vice versa, where 

( 6 - 7 ) S ■= [_b o) ' S ~ l = \B~ 1 -B-^AB- 1 

Proof. Property (J6.6)) follows immediately from the relation below, which we obtain from 
integration by parts: 

W*SW\H — {{L — \y(j, U) L i {xuX2) - (U, (L - X)U) LHxuX2) = 0, 

by the definition of the adjoint operator. (Indeed, the righthand side may be viewed as 
defining the quadratic form S.) □ 

6.2. Domain of consistent splitting. Define 

(6.8) A:=nA±, j = l,...,n, 

where denote the open sets bounded on the left by the algebraic curves A^ (£ ) determined 
by the eigenvalues of the symbols —£, 2 B± — z£A± + C± of the limiting constant-coefficient 
operators 

(6.9) L±U := B±U" - A±U' + C±U 

as £ is varied along the real axis. The curves A^(-) comprise the essential spectrum of 
operators L±. For Majda's model, the computations of Section l4~Tl yield 

(6.10) A C Q v := {A : KA > max{-r/i|3A|, r/ 2 |3A| 2 }}, rj, > 0. 

Lemma 6.2 QMaZ3 ). The set A is equal to the component containing real +oo of the domain 
of consistent splitting ( defined in Section \4-l\j for the eigenvalue equation of L written as a 
first-order ODE (pTT|) . 

6.3. Basic solution formula. Let 

(6.11) $+:=P + tf, J = l...,n, 
and 

(6.12) $t : = p_yr, j = n + 1, . . . , 2n 

denote the locally analytic bases of the stable manifold at +oo and the unstable manifold at 
— oo of solutions of the eigenvalue equation ()4.6|) written as a first-order system (|5.1|) that 
were found in Section I5~2l (i.e.. the normal modes), and set 

(6-13) *+:=($+...,$+), $-:=($; +1 ,...,$ 2 - n ), 

and 

(6.14) $:=($+,$-). 

Define the solution operator from y to x of ()4.6|) , denoted by jF y ^ x , as 

(6.15) T y ^ x = <t>{x,X)<t>- 1 {y,X) 
and the projections ILj 1 on the stable manifolds at ±oo as 

(6.16) II+ = ($+(j/,A) 0)<S>-\y,\) and U~ = ( ®~(y, A) ) A). 
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Then, we have the following general result established in [ZH, MaZ3 . 

Proposition 6.3 ffZHl IMaZHj ). With respect to any L p , 1 < p < oo, A consists entirely 
of normal points of L, i.e., resolvent points, or isolated eigenvalues of constant multiplicity. 
On this domain, the resolvent kernel G x is meromorphic, with representation 



(6.17) G x {x,y) 



{I n ^)^ x Ii+S-\y)(I n ^) tr x>y, 
- ( J n , 0)^* n- S-\y) (I n ,0) tr x<y, 

S^ 1 as described in (16. 7|) . Moreover, on any compact subset K of p(L) fl A (p(L) denoting 
resolvent set), there hold the uniform decay estimates 

(6.18) \dld k y G x (x,y)\<Ce~^-y\, 

< \j\ + | A: | < 1, where C > and f] > depend only on K , L. 

Remark 6.4. Formula (j6.17|) extends |ZH|. IMaZ3] to the full phase-variable representation 



(6.19) 



G x (x,y) {d/dy)G x {x,y) \ /jR-n+S"^) x > y, 

(d/dx)G x (x,y) (d/dx)(d/dy)G x (x,y)J U~ S^y) x<y. 



6.4. Generalized spectral decomposition. Formula (|6.17|) suffices for the description 
of intermediate- and high-frequency behavior. For the treatment of the key low-frequency 
regime, it is preferable to use a modified representation of the resolvent kernel consisting of 
a scattering decomposition in solutions of the forward and adjoint eigenvalue equations. 

From ()6.6|) . it follows that if there are n independent solutions . . . , 0+ of (L—XI)U = 
decaying at +oo, and n independent solutions . . . , <t>2n °f ^ ne same equations decaying 
at — oo, then there exist n independent solutions ipn+ii • • • > 4>tn °f (L* — X*I)U = decaying 
at +00, and n independent solutions ipi , . . . , ip~ decaying at —00. Precisely, setting 

(6.20) #+:=P + V/, j = n + l...,2n, 

(6.21) *J:=P_Vf, j = l,...,n, 
and 

(6.22) 

similarly as in (|6.11|) - (|6.14p . where ^/f = (ipf, (tpf)') are exponentially growing solutions ob- 
tained through Lemma IS~T1 we may define dual exponentially decaying and growing solutions 
tpf and (f>f via 

(6.23) ®)± = I- 
Corollary 6.5 ( |ZTllMaZ3j ). On An p(L), there hold 

(6.24) G x {x,y) = ^ M+(A)0+(x; A)fe(y; A)* 

/or y < < x, 

(6.25) G A (x, y) = J2 ^ k (X)<l>j(x; A)fe(y; A)* - A)^~(y; A)* 
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for y < x < 0, and 



(6.26) 



G\(x,y) 



d jkW<Pj& (y; a)* + x )tk(y> A ) 



j,k k 



for x < y < 0, with 



(6.27) 



M 



(-1,0) ($+ $-) * 



and 




Symmetric representations hold for y > 0. 



Proof. Rearrangement of (J6.17)) using (J6.23j) ; see [MaZ3j. 



□ 



Remarks 6.6. 1. This representation reflects the classical duality principle (see, e.g. [ZH], 
Lemma 4.2) that the transposition G* x (y,x) of the Green's function G\(x,y) associated with 
operator (L — A) should be the Green's function for the adjoint operator (L* — A*). 

2. In the constant-coefficient case, with a choice of common bases \P = $ T at ±oo, 
()6.24|) - ()6.28|) reduce to the simple formula 



This reveals an analogy to the usual representation obtained by Fourier transform solution. 
We see that ()6.25|) (resp. (j6.26|0 in the far- field limit x — > +oo (resp. x — > — oo) consists of 
the limiting constant-coefficient resolvent kernel plus an exponentially decaying error term. 

6.5. Resolvent kernel bounds. The basic bound (j6.18|) is sufficient to treat intermediate 
frequencies r < \X\ < R, r, R > (indeed, it is difficult to say more in this regime). As 
|A| — ► oo, orA^Oe dp(L), however, the bounds are not uniform, and so separate analyses 
are needed in these high- and low-frequency regimes. 

In the high-frequency (~ short-time) regime, we have the following classical analytic 
semigroup-type bounds following from strict parabolicity alone. These may be obtained 
from the basic formula (J6.17j) using the parabolic rescaling x — > a;|A| 1//2 and estimates (the 
Tracking Lemma of |ZHIZ2trMaZ3p for slowly-varying-coefficient ODE; see, e.g., |Sa|IZHtlZT] . 

Proposition 6.7 (High-frequency bounds [ZHJ). For general strictly parabolic systems of 
form (|6.3|) . R, C > sufficiently large, and rji, 7] 2 , 9 > sufficiently small, 



(6-31) \G x (x,y)\ < C \\\-^e- e ^-y\, \d x G x (x,y)\, \d y G x (x,y)\ < Ce^ x ^ 



for all A £ £l v \ B(0, R), with Q v as defined in (j6.1()j) . 

Thus, the only resolvent bounds that depend on the details of the model are the crucial 
low-frequency (~ large-time) bounds carrying information relevant to large-time asymptotics. 
For the Majda model, these are as follows. 
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(6.29) 




x > y, 
x < y, 



where, generically, <f> • , <f> • may be taken as pure exponentials 



(6.30) <ft{x)$?{v) = e^ (A)(:c - y) Vf (A)Vf *(A). 



Proposition 6.8 (Low-frequency bounds). Let U = (u,z) be a traveling wave profile of 
Majda's model, satisfying (D). Then, for r > sufficiently small, the resolvent kernel G x 
has a meromorphic extension onto 5(0, r) C C, which may be decomposed as 

(6.32) Gx = E x + S x + R x , 
where 

(6.33) E x (x,y) : (a-;^^,-^^^— 2 ^ 3 ^ «" > 0, 
V ' K ' \\- 1 U'(x)7r-(y) tI or < 0, 

for y < 0, with irj bounded solutions of the adjoint eigenvalue equation for X = 0, nj 
convergent as y — > — oo to cLj and 7r f 7 exponentially decaying as y — > — oo, 



(6.34) E x (x,y):={ 



r 1 U'(x)7r+(y) tI e^ a+ - x2 / a+3 ^ 
+ \- 1 U'(x)iT+(y) tT e ( -- x / s+x2d / s3 ^ 



\- 1 U'(x)7r+(y) tT e^ s+x2d / s3 ^ a+ > 0, 



fory > 0, with irt bounded solutions of the adjoint eigenvalue equation for X = 0, convergent 
as y — > +oo; 



(6.35) S x (x,y) :-- 



Cj R+L7*e(- A /- + +^ 2 /« +3 )-+(A/a--A 2 /«- 3 ) S / a - a + > o 



otherwise 
/or y < < x, Rj and Lj constant vectors as defined in (|4.24j) . 



(6.36) S x (x,y) :-- 

for y < x < 0, and 



: tr e (_A/a-+A 2 /a- 3 ) 

otherwise 



R - L f e i-^h~+*/^){z-y) a - > o 



(6-37) 5 A (a;,y):=^ } 1 . 

I otherwise 

for x < y < 0, with similar relations for y > 0; and R x denotes a faster- decaying residual 
term. 4 

Remarks 6.9. 1. Recall, (ET2U1) - @2Z) , (jOifl - ljOHf . that vectors = (*,0) tr appearing 
in scattering terms S\ have vanishing ^-component, and also Lp _L (— g, 1), a fact that will 
be important in our later nonlinear stability analysis. 

2. The case a~ < in (j6.33J) occurs only in the extreme situation of a strong deflagration, 
for which there are no incoming characteristics on the lefthand side y < 0. This is essentially 
the only difference from the corresponding proposition for viscous shock waves in the general 
Lax or undercompressive case, and represents just an anomaly in bookkeeping. 



See, e.g., |MaZ3l 172] for bounds in the viscous Lax shock ~ strong detonation case. 
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Proof. Similarly as in the viscous shock case treated in [MaZ3j, this follows from represen- 
tations ()6.24j) - (j6.28j) of Corollary 16.51 estimating modes <fif, ipj , (f)f, and ipj using the as- 
ymptotic description given by the Conjugation Lemma together with the constant-coefficient 
analysis of Section 14.11 and estimating scattering coefficients Mj kl d- k crudely by Laurent 
series: e.g., 

djk = dj^X 1 + dj k = d® k + . . . , 

noting that pole terms of order A~ fc correspond to zeroes of order k of the Evans function, 
hence (by (£>)) are at most order k = 1 and (without loss of generality coordinatizing so that 
'Pi = ^2n = U\x)) involve only zero-eigenfunction U'(x) as x-dependent factor. Specifically, 
E\ comprises exact pole terms, while S x comprises order one terms involving products of 
slowly decaying forward and dual modes (i.e., modes that are merely bounded for A = 0), 
the latter estimated to exponentially decaying error via the Conjugation Lemma, while R x 
comprises remaining, residual terms. 

Vectors Rj, Lj in the formulae for S x arise through the limiting, constant-coefficient anal- 
ysis of Section 14.11 Finally, the information that ttJ at A = converges as y — ► — oo to 
cL~f if ct > and to zero if a" follows by the fact that in the first case there exists but a 
single bounded, nondecaying solution of the eigenvalue equation as y — ► — oo, with asymp- 
totic direction Lj, and in the second case there exists no bounded, nondecaying solution. 
These facts, in turn, are readily verified by the constant-coefficient analysis of Section 14. 1[ 
combined with the Conjugation Lemma. □ 



6.6. The system case. Proposition 16.51 admits a straightforward generalization to the sys- 
tem case (i.e., the vectorial version of (|1.3|) discussed in Remark |2~T1 3). substituting in place 
of a ± 1 Ly, R^ the eigenvalues aj and eigenvectors Lijp ■ of f u (U±) — si and in place 
of L+, Rf the eigenvectors L^, R^ { of the (now matrix- valued) diffusion coefficient d, as 
described in Section 14.21 (recall, there are no slow reactive modes on the minus infinity side 
x < 0). See Proposition 4.22, |Z2j for a corresponding description in the viscous shock case. 

Proposition 6.10 (Low-frequency bounds). Let U = (u,z) be a traveling wave profile of 
vectorial Majda's model under dissipativity hypothesis (J3.15j) . satisfying (D). Then, forr > 
sufficiently small, the resolvent kernel G\ has a meromorphic extension onto B(0,r) C C, 
which may be decomposed as 

(6.38) G X = E X + S X + R x , 

where 



(6.39) E x {x,y):-- 



\- l U' {x)ti; {y) tv alla7<0 



for y < 0, with n bounded solutions of the adjoint eigenvalue equation for A = ; expo- 
nentially convergent as y — > — oo to cj ■ and ir~ i exponentially decaying as y — > — oo ; 



'Recall, Remark H1.7|l . this implies in part that U is a unique, transversal connection. 
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aj<0 



(6.40) E x (x,y):={ 



some ol- < 0, 



{-X/s+\ 2 d+/s 3 )y 



i=l 



A" 1 J2 U\x)n+{y) tr e { - X/s+X2dt/s3)y all a+ > 0, 



/or y > ; with tt^j, 7r+ bounded solutions of the adjoint eigenvalue equation for A = ; 
exponentially convergent as y — > +oo; 



(6.41) S x (x,y):= £ c^J^L^"*/" 



, +A 2 b+/a+>+(A/aT-A 2 6-/a7°) 2 , 



a fc ,aT>0 



/or y < < x ; i?^^ and Lj fc constant vectors as defined in Section 



(6.42) 



/ j Ji3 Jft 



a f,3< > a k>° 



for y < x < 0, and 



S\(x,y) := R k L k e 



(6.43) 



-r -* J-X/a^+X%/a^)(x-y) 



«7,i <0 ' Q fe >0 



/or x < y < ; where , c^' ± are scalar constants, with similar relations for y > 0; ana 1 i?A 
denotes a faster- decaying residual term. 



7. Green function bounds 

We may now estimate the Green function G(x,t;y) := e Lt d~ y (x) associated with the lin- 
earized operator L about the wave, determined by 

(7.1) (d t -L)G = 0, G(x,0;y) = 5 y (x), 

via the inverse Laplace-transform formula, following the approach of [ZH, lMaZ3] . We present 
our results using a bookkeeping scheme similar to that of HZJ in the undercompressive 
viscous shock case. 
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7.1. Basic bounds. Recall the standard notation 



erring) : = 




Proposition 7.1. Let U = (u, z) be a traveling wave profile of Majda's model, satisfying 
(D). Then, the Green function G(x,t;y) associated with the linearized equations ()6.3|) may 
be decomposed as G = E + G, where 

(7.2) E(x,t;y) = U'(x)e(y,t), 



(7.3) e(y,t) 



for y < and 



(7.4) e(y,t):={ 



errfn ' v+a 1 



v r (y) 



errfn 



y + a + t 



^fn ( ltjpJJ ^(y) a > 0, 

a" < 0, 



At 



errfn 



y + a + t 



r I y + st\ 
errfn — — errfn 



'Adt J 



At 
y- st 

VAdt 



T?(y) 



K(y) 



k (err/n - err/n (^=§)) 7T+(y) 

/or > ; with nf as in Proposition ^.^ in particular, 



a< < 0, 



a+ > 0, 



\*f\<c, 



\d v nf\ < Ce-*b 



(7.5) 

with \ttJ \ < Ce~ v ^ if a~ < (strong deflagration case) andn^ — > as x — * — oo otherwise, 

and, denoting by af the "undamped" characteristic speeds aj = a~ , = a + , and a+ = — s 
(Note: a~ does not appear), 

\d« y G(x,t;y)\<Ce-^-y^ 

+ C(r H/2 + \a x \e~^ + \a y \e- vlyl ) 



x | t -l/2 e -(x~y~a-t) 2 /Mt e - V x+ 
k 



(7.6) 



a fc >0, o 3 . <0 



+ Yl X{K*[>ll/[}* 



(x-a+(t-|y/a-|)) 2 /Mt^~^- 



a k >0,aj>0 



(i+ S (t-|y/a-|)) 2 /Mt -rjlccl 



± 



< I or | < 1 for y < and symmetrically for y > 0, for some rj, C , M > 0, where x 
denotes the positive/negative part of x and indicator function X{\ a ~t\>\ y \} ^ s 1 f or \ a kA — \v\ 
and otherwise. Moreover, for x < 0, |(0, l)G(x, t; y)\ decays at the faster x-derivative rate 
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a x = 1, as does | (0, l)G(x, t; y)(l, 0) tr | for any x, and, for y < 0, \G(x, t; y){— q, l) tr | decays 
at the faster y-derivative rate a y = 1. 

Proof. Reflecting the formal relation that G\ is Laplace transform of G, we have the Inverse 
Laplace transform formula 



(7.7) G(x, t;y) = ^~ f e M G x (x, y) dX, 

where T = d{\ : ReX > Q\ — #2|^A|} is the boundary of an appropriate sector containing 
the spectrum of L, 6 2 > 0. Following |ZHj IMaZ3j . we may thus convert the detailed resolvent 
kernel estimates of Proposition 16. 51 to estimates on the Green function via stationary phase, 
or Riemann saddlepoint, estimates on (|7.7|h exactly as was done in the viscous shock case. 

Specifically, using the property that G\ is meromorphic on Q v U 5(0, r) for r, 77 > 
sufficiently small (see Propositions 16. ^ 16. 7| and Ifi.Hj) and analytic on the the resolvent set 
p(L), we may estimate the contribution of each of the various meromorphic components 
C\ of G\ by a combination of direct evaluation using Calculus of residues and strategic 
deformation of the contour T so as to minimize 

\tn „,M \a\ \ — 4, ^KAt 



e M C x (x,y)\\d\\= d> e uxt \C x {x, y)\ \dX\ 
v Jr 

for each fixed x, y, t, with the main contribution to E coming from explicit evaluation of 

the corresponding low- frequency term E\ in Propositon 16.81 and the main contribution to G 

coming from explicit evalution of S\. See [ZHl IMaZ3l IZ2j for details. 

It remains only to verify the key properties of faster decay of |(0, l)G(x, t; y)\ for x < 0, 

|(0, l)G(x, t; y)(l, 0) tr | for general x, and \G(x, t; y)(— q, l) tr | for y < 0. To see the first 

property, we have only to observe, in the bounds of Proposition 16. 5[ that, for x < 0, only 

fluid modes appear in the rate-determining term S\, and these lie in direction Rj = (*, 0) tr 

having vanishing z-component. 6 For x > 0, fluid modes again lie in direction Rlf = (*,0) tr 

having vanishing ^-component, while reactive terms appear as scalar multiples of projector 

R+(L+) tr = (0, *) orthogonal to (1, 0) tr . Thus, |(0, l)G(x, t; y)(l, 0) tr | to lowest order involves 

only ^-components of fluid terms, hence again is faster decaying; this yields the second 

property. Likewise, for y < 0, S\(—q, l) tr = 0, since _L (—q, l) tr , and this yields the third 

property, completing the proof. □ 

Remarks 7.2. 1. Similarly as in the viscous or relaxation shock case, the bounds of 
Proposition 17.11 may be interpreted as describing the evlution of an initial delta-function 
perturbation at y as the superposition of signals convecting along hyperbolic characteristics 
and diffusing as approximate Gaussians until they strike the shock layer, whereupon they 
scatter as reflected and transmitted waves along outgoing characteristics, at the same time 
exciting the stationary mode U' . The main new feature in the reacting as compared to the 
nonreacting case is the exiting signal along the reaction characteristic on the lefthand (x < 
0) side, for which a constant-coefficient analysis indicates that the Gaussian signal is now 
exponentially decaying in time, due to burning of the reactant. This is reflected in the final 
term of f)7.6j) . consisting of an ordinary Gaussian reflected left into an exponentially penalized 
field e~ 9 ^, a term indistinguishable in modulus bound from a Gaussian multiplied by a factor 



6 The restriction x < is necessary because of incoming (i.e., leftmoving) undamped reaction waves for 
the case y > not listed in Proposition 17. II 
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decaying exponentially in the travel time after reflection. We call this leftgoing reactive 
characteristic speed "damped" and all others "undamped" to distinguish this behavior. 

2. A second difference between the reacting and nonreacting case, this time confined to 
(scalar) Majda's model, is the different structure of the excited term E in the case a' < 0, 
a book-keeping anomaly arising because of the absence of incoming waves on the lefthand 
side in this (strong deflagration) case. This different structure has essentially no effect on 
the analysis; see Remark 18.51 

3. The improved bounds for |(0, l)G(x, t; y)\ and \G(x, t; y)(—q, l) tr | for x < are similar 
to those of the relaxation case [MaZlJ, to which the linearized equations are analogous on 
the side x < on which <p > 0, and play a similarly important role in the later nonlinear 
stability analysis. In conservative coordinates w := u + qz, z of (|2.2|) - (|2.3j) . these bounds 
have the simpler statement that the Green function decays more rapidly in its z-components, 
both output and input. 

7.2. The system case. The somewhat cumbersome summation notation of Proposition [7^ 
is designed for easy generalization to the system case. Indeed, starting from Remark 16.91 3. 
it is straightforward to verify the analogous theorem for the full reactive Navier-Stokes 
equations with artificial viscosity- more generally, the abstract vectorial model u G M. n , z 6 
]R m described in Remark 12.11 3- with undamped characteristic modes now af = af, . . . , a^, 
—s (— s multiplicity m = dimz), af and Lrf denoting the eigenvalues and left eigenvectors 
of (df /du){U±) — si. See [HZJ for the analogous description in the viscous shock case. 

Proposition 7.3. Let U = (u,z) be a traveling wave profile of vectorial Majda's model 
under dissipativity hypothesis (J3.15|) . satisfying (T>). 7 Then, the Green function G(x,t;y) 
associated with the linearized equations ()6.3|) may be decomposed as G = E + G, where 

(7.8) E(x,t;y) = U'(x)e(y,t), 



(7.9) e(y,t)= l Ea I>° 

for y < and 
(7.10) 



errfn 1 l 



Ubjt 



errfn 1 



46" t 



e{y,t) 



j errfi 



11 



y + ajt 



Abjt 



errfn 



y + ajt 



+ ( errf 



y + st \ ( y-st 

n | — , I — errfn 



8=1 



7Tf j(y) tr some a- > 0, 



all a ■ < 0, 



some at < 0, 



tr 



Ei=i (err/n (^=T^J - err f n ( 



y-st 



all aj > 0, 



Recall, (jl.7|l . this implies in part that U is a unique, transversal connection. 
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for y > 0, with as in Proposition ^.^ in particular, 
(7.11) \*?\<C, \d v *f\ < Ce-*M, 

with \ < Ce~ v ^ ifa~ < (strong deflagration case) andirjr — > as x — > — oo otherwise, 
and G satisfies ([7.6|) . Moreover, forx < 0, |(0, I r )G(x, t; y)\ decays at the faster x-derivative 
rate a x = 1, as does |(0, I r )G(x,t;y)(I n ,0) tT \ for any x, and,fory<0, \G(x,t;y)(—q tv ,I r ) tr \ 
decays at the faster y -derivative rate a y = 1. 

7.3. Linearized stability criterion. 

Proof of Theorem \l.°A Sufficiency of {T>) for linearized orbital stability follows immediately 
by the bounds of Theorem 17.11 (resp. Remark 17.21 4) and standard L q — > L p convolution 
bounds, exactly as in the viscous shock case, setting 

/+oo 
E(x,t;y)u (y)dy 
-oo 

so that 

/+oo 
G(x,t;y)u (y)dy; 
-oo 

sec [ZH, MaZ3j IZ2j for further details. Necessity follows from more general spectral con- 
siderations not requiring the detailed bounds of Theorem 17-H see the discussion of effective 
spectrum in |ZH| IMaZ3| IZ2j . The argument goes again exactly as in the viscous shock 
case. □ 



8. Nonlinear stability 

We can now readily establish nonlinear stability by a combination of the methods used 
in |HZj to treat general undercompressive viscous shock waves (~ x > behavior), and 
the methods used in |MaZlj to treat relaxation shocks (~ x < behavior). As it costs no 
additional effort in bookkeeping, we carry out this part of the argument in the full generality 
of the system case. Recall fRemark l2.H 3) that this includes the artificial viscosity version of 
the full reactive Navier-Stokes equations with multi-species reaction and reaction-dependent 
equation of state. 

Denoting by af the "undamped" characteristic speeds: aj = a~ , = a + , and a+ = — s 
for Majda's model; aj ri = a~ , i = 1, . . . , n, i = af , i — 1, ...,n, = — s, i — l,...,min 
the system case u G M™, z G M m ; define 

(8.1) 9(x, t) := ^ (1 + t )-i/2 e -l— 7*l 2 /^ + (1 + t)- 1 ' 2 e-\ x - a M 2/Lt , 

aj<0 a+>0 

iptix, t) := x (x, t) C 1 + \ x \ + tY ll2 (X + \x- ajt\)- 1/2 

aj<0 

+ X(x, t) J2 I 1 + \ x \ + ty 1/2 (l + \x- ap\y 1/2 , 

a+>0 
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(8.2) 



and 

M*> *) ■= (! " X(x, t))(l + \x- a-t\ + *V2)-3/2 
+ (l- X (x,t))(l + |x-a+t|+t 1 /2 ) -3/ 2) 

where L > is a sufficiently large constant and £) = 1 for 

(8.4) x G [min{a~t, 0}, maxja^t, 0}], 

3 3 

that is, for x between the extremal outgoing undamped characteristics, and zero otherwise. 8 
Then, we have the following pointwise version of Theorem 11.51 

Proposition 8.1. Let U(x — st) be a traveling combustion wave of Majda's model (more 
generally, the vectorial version including reactive Navier-Stokes equations with artificial vis- 
cosity) and \Uq(x)\ < E (l + \x\)~ 3 ^ 2 , E sufficiently small. Then, there exist S(-) and 5(+oo) 
such that 

\U(x,t) - U m (x)\ < CE (6 + + ih)(x,t), 

(8.5) \5(t)\<CE (l + t)-\ 

\5(t) - 5(+oo)\ < CEoil + t)- 1 / 2 , 

where U denotes the solution of the same equations with perturbed initial data Uq — U + Uq. 

As discussed in the introduction, we establish Proposition 18.11 by a combination of the 
analysis of undercompressive viscous shock waves in |HZj and of relaxation shocks in [MaZlj. 
Following |HZj . set 

(8.6) U(x,t):=U(x + 6(t),t)-U(x), 
so that (|6.1|) becomes by Taylor expansion of F, G: 

(8.7) U t -LU = Q(U) X + R(U) + 8(t)(U x + U x ), 
L as in (|6.3|) . where R{U) = (—q, l) tT r(U), r(U) scalar, with 

r(U) = 0{\U\ 2 e- rfB ), 
so long as \U\ remains bounded, where x + denotes the positive part of x and rj > 0. 
Remark 8.2. Here, in the description of R, r we have used the specific form 

(8.9) G(U) = -0(«) H k 
of the reactive source in (|6.1|) . together with Taylor expansion 

(8.10) (<j)(u + u)(z + z) - (4>(u)(z) - {(j)'{u)uz + (f)(u)z) = (p'(u)uz + <f>"(u + 6u)u 2 z, 

< j3 < 1, and the fact that 4>'{u + v) < Ce~ vx+ for |u| sufficiently small, by assumption 
u + £ [ui,u l ], the property that 4>'{u) = for u G" and exponential convergence of 

U(x) to U+ as x — ► +oo. This computation, and its exploitation in the later argument (see 

8 This repairs a minor omission in |HZ| . where H8.4[) was stated incorrectly as x € [miiij{a~t}, maxj{a^"i}]. 
The formulae differ in the case that there are no outgoing characteristics on one side: extreme Lax shock or 
strong detonation. 
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especially the auxiliary bounds of Lemma [8 .6J) . are the main new features in the combustion 
context as compared to the undercompressive viscous shock wave case. 

Recalling the standard fact that U' is a stationary solution of the linearized equations 
(P , LtJ' = 0, or 



G(x,t;y)U x (y)dy = e Lt U x (x) = U'(x), 
we have by Duhamel's principle: 

POO 

U(x,t) = I G(x,t;y)U (y)dy 

G(x, t - s; y)(-q, lfr(U)(y, s) dyds 
G y (x,t- s; y)(Q(U) + 5U)(y, s) dyds + 8(t)U'(x). 

JO J~oo 

Defining 



J — oo 

t POO 



</-oo 
t poo 



5{t) = - e(y,t)U (y)dy 



oo 

t POO 



8.11) -/ / e{y,t-s){-q,l) tr r{U){y,s)dyds 



J-oo 

t P + OO 



+ / / e y {y,t-s){Q{U)+6U)(y,s)dyds, 



y 

J-oo 



following fZal IMaZll IMa72l IMaZ4j . where e is defined as in CTSD-dHJ) (that is, e = £\ e s ) 
and recalling the decomposition G = E + G, we obtain finally the reduced equations: 



U(x,t) = G(x,t;y)U (y)dy 

J — OO 

rt roc 

(8.12) +/ / G(x,t-s;y)(-q,lfr(U)(y,s)dyds 

JO J-oo 



t poo 

G y (x,t - s;y)(Q(U) + 5U)(y, s)dyds, 



J-oo 



and, differentiating (|8.11|) with respect to t, and observing that e y (y,s) —r as s — >• 0, as 
the difference of approaching heat kernels: 



5(t) = - e t (y,t)U (y)dy 

J —oo 
pt poo 

(8.13) + / / e t (y,t-s)(-q,lfr(U)(y,s)dyds 

JO J-oo 
pt P+OO 

+ / / e yt (y,t-s)(Q(U)+SU)(y,s)dyds. 

JO J-oo 



The following integral estimates are established in HZJ. 

32 



Lemma 8.3 (Linear estimates |HZj ). Under the assumptions of Theorem \1.5l 



oo 

+00 



(8.14) 



\G(x, t; y)\(l+ \y\)- 3 / 2 dy < C{6 + fa + fa) {x, t), 
e t (y,t)\(l + \y\)- 3 / 2 dy<C(l + t)-^, 
(y,t)\(l + \y\r 3 ^dy<C, 
t(U.t)- e(y, +00) I (1 + M)- 3 / 2 dy < C(l + t)^ 2 , 



OO 

+00 

\e[ 

00 



for < t < +00, some C > 0, where G and e are defined as in Proposition 7.1 



Lemma 8.4 (Nonlinear estimates HZJ). Under the assumptions of Theorem ] 1.51 



t i>+oo 

/ \G y (x, t - s; y)\V(y, s) dyds < C(6 + fa + ^2) fa t), 

J -00 

/ / \e yt (y,t-s)\q(y,s)dyds<C(l+t)- 1 , 

J0 J -00 



[8.15) 



00 r+00 



\e y {y, +oo)|*(y, s) dy < C T (1 + t)" 1 / 2 , 



—00 



n f+00 

/ / \e y (y,t-s)-e y (y,+oo)\y(y,s)dyds<C(l+t)- 1 / 2 , 

J0 J -oo 

for < t < +00 , some C > 0, where G and e are defined as in Proposition ]! '. 1\ and 

*(!/, s) := (1 + sY^s-^iO + fa + faf(y, s) 

{Kb +(i + s)- 1 (0 + fa + fa)(y,s). 

Remark 8.5. The case a(a~) < in (|7.3|) . (|7.9p . occurring for strong deflagrations as 
described in Remark 16.91 2. is the only one requiring discussion, since in all other cases 
the bounds are identical to those of the shock case. We have only to note that |e|, \e y \, 
\zyt\ = =0 in this case also satisfies the same bounds (or better) that are actually used 
in the proofs of Lemmas 18.31 and 18.41 

To these, we add the following auxiliary estimates special to the combustion case. 
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Lemma 8.6 (Auxiliary estimates). Under the assumptions of Theorem \l.b\ 
(8.17) 

"t f+OO 

\G(x, t - s; y){-q, lfe~™ + \V(y, s) dyds < C(9 + fa + fo)(x, t), 

'0 J -oo 

ft r+oo 



rt r>+oo 

/ / \e t (y,t-s)(-q,l) tI e-™ + \*(y,s)dyds<C(l+t)- 1 

JO J-oc 

/■+oo r+oo 

/ / \e(y, +oo)(-q, lfe~™ + \*(y, s) dy < C7 7 (l + t)" 

</-oo 

/ / I (e(y, t - s) - e(y, +oo)) (-g, l) tr e"™ + |*(y, s) dyds < C(l + t)' 1 

JO J-oc 



1/2 



for < t < +oo ; some C > 0, G and e as in Proposition \ 7. 1\ and \I/ as in f!8.16|) . 

Proof. For the G-estimate, we have only to recall that, by the bounds of Proposition 17.11 
G(—q, l) tr obeys the bounds of G y for y < 0, while Ge~ m+ obeys the bounds of for 
y > 0. The e-estimates follow similarly, by the observation that, for y < 0, e(—q, l) tr and 
e t (—q, l) tr decay like (e~ r? ' y ' + (l+t)~ 1 / 2 ) times the bounds for e and e t , hence obey the bounds 
for \e y \ and e yt , since tcJ is asymptotically parallel to ± (— g, l) tr , with convergence at 

exponential rate. Likewise, for y > 0, ee~ m+ and ete~ vy+ decay like (e _,? ' y ' + (1 + t) -1 / 2 ) 
times the bounds for e and e t , hence obey the bounds for \e y \ anb \e yt . Thus, all bounds 
follow by the same arguments as in the proof of Lemma 18.41 □ 



Proof of Proposition \8.1l With these observations, the proof of nonlinear stability goes es- 
sentially as in [HZ . Define 

(8.18) C(*):= sup (V|(0 + Vi +ik)-\y,t) + |<j(s)|(l + s)Y 

y,0<s<t \ ' 

We shall establish: 

Claim. For all t > for which a solution exists with ( uniformly bounded by some fixed, 
sufficiently small constant, there holds 

(8.19) at)<C 2 (E + ((t) 2 ). 



From this result, provided E < l/4Cf, we have that ((t) < 2C 2 E Q implies ((t) < 2C 2 -E'cb 
and so we may conclude by continuous induction that 

(8.20) C(*) < 2C 2 E 

for all t > 0. (By standard short-time existence for artificial viscosity systems (see, e.g., 
[HoSl IZTI] ). U G C 1 exists and ( remains continuous so long as ( remains bounded by some 
uniform constant, hence ()8.20|) is an open condition.) Thus, it remains only to establish the 
claim above. 

Proof of Claim. We must show that U{6 + ipi + V^) -1 and |<j(s)|(l + s) are each bounded 
by C(E + ((t) 2 ), for some C > 0, all < s < t, so long as £ remains sufficiently small. 

34 



Recalling definition (|8.18|) . we obtain for allt > and some C > that 
(g21) \5(t)\<((t)(l + t)-\ 

\u{x,t)\ <c(t)(0 + Vi + ^)(^*), 

and therefore 

(8.22) \(Q(U) + 5U)(y,s)\ < C((t) 2 *(y,s) 

with \1/ as defined in (J8.16|) . for < s < t. 

Combining (J8.22JI with representations (|8.12|) - (j8.13|) and applying Lemmas 18.31 and 18. 4[ 
we obtain 

/oo 
\G(x,t;y)\\U (y)\dy 
-oo 

/»£ poo 

+ / / |e(y,*-s)(-g,l) te ||r(l7)(y,s)|dytfe 

J </-oo 
/»£ /*oo 

+ / / |G,(x,t- S ;y)||(g(i7) + ^)(y, S )|dyd S 

JO ./-oo 



/oo 
|G(x,t;y)|(l + |y|)- 3 / 2 ^ 
-OO 

/*£ /*oo 

CC(t) 2 / | e (y,t- a )(- gi l) te e-^|*(y,a)|c^cia 

JO ./-oo 



'0 J-oo 

ft poo 

\2 



+ C((t) 2 / |G^(z,*-s;y)|tt(j/,s)dj/ds 

J0 ./-OO 

<C(E + C(t) 2 )(6 + ^+?p 2 )(x,t) 



and, similarly, 

< / \e t (y,t)\\Ua(y)\dy 



OO 

rt poo 

+ / / \e t (y,t-s){-q,ir\\r{U)(y,s)\dyds 

JO J-oo 

ft P+OO 

+ / / \e yt (y,t-s)\\(Q(U)+5U)(y,s)\dyds 

Jo J-oo 

oo 



< / E \e t (y,tm + \y\r 3/2 dy 



-oo 

ft p+oo 

+ / / CC(t) 2 \e t (y,t-s)(-q,ire-™ + \*(y,s)dyds 

Jo J-oo 



ft p+oo 

/ / Ct(t) 2 \e vt (y,t-s)\V(y,s)dyds 

Jo J-oo 



+ 



<c(£ +cw 2 )(i+tr :l 

Dividing by (8 + ipi + ip2){%-i t) and (1 + £) _1 , respectively, we obtain ()8.19j) as claimed. 
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From (|8.19jh we obtain global existence, with ((t) < 2CEq. From the latter bound and 
the definition of £ in ()8.18|) we obtain the first two bounds of (|8.5|) . It remains to establish 
the third bound, expressing convergence of phase S to a limiting value 5(+oo). 

By Lemmas I8.3H8.4I together with the previously obtained bounds (|8.22|) and ( < CE , 
and the definition (J8.18)) of (, the formal limit 

/oo 
e(y,+oo)U (y)dy 
-OO 

rt poo 

+ / / \e(y,t-s)(-q,l) tr \\r(U)(y,s)\dyds 

JO J-oo 
r+oo r+oo 

+ / / e y (y,+oo)(Q(U) + 5U)(y,s)dyds 

JO J-oo 

/oo 
\e(y,+oo)\E (l + \y\)- 3 / 2 dy 
■oo 
rt r+oo 

+ / \e(y,t-s)(-q,l) tt e-™ + \CE *(y,s)dyds 

Jo J-oo 



-oo r+OO 

+ / / \e y (y,+oo)\CE ^f(y,s)dyds 



< CE, 



o 



is well-defined, as the sum of absolutely convergent integrals. 
Applying Lemmas I8.HH8.4I a final time, we obtain 

/oo 
\e(y,t) - e(y,+oo)\\U (y)\dy 
-oo 

rt poo 

+ 



J-oo 
t r+oo 



J-oo 

OO r+OO 



e(y, t- s) - e(y, +oo) J (-q, l) tr | \r(U)(y, s)\dyds 
\e y (y,t- s) - e y (y, +oa)\\(Q(U) +5U)(y,s)\dyds 
\e y (y, +oo)\\(Q(U) + 5U)(y, s)\dyds 
< I \e(y,t)-e(y,+oo)\E Q (l + \y\)- 3 / 2 dy 

—oo 

+ lo r 1 ( e(y ' t ~ s) ~ e(y ' +oo) ) 1 ) tre ~ nu+ \ CE ^ s ) d y ds 

rt r+oo 

+ \e y (y,t- s) - e y (y,+oa)\CE ^(y,s)dyds 

Jo J-oo 
r+oo r+OO 

+ / / \e y {y,+oo)\CE Q V(y,s)dyds 

Jt J-oo 



<CE (l + t) 



-oo 
-1/2 



establishing the remaining bound and completing the proof. □ 
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Proof of Theorem \l.b\ Immediate from Proposition 18.11 by integration of bounds ()8.5|) . □ 

Remark 8.7. Proposition 18.11 gives a time-asymptotic description of perturbation U as a 
superposition of algebraically decaying signals propagating along outgoing undamped char- 
acteristic directions. A brief examination reveals that these consist entirely of fluid dynam- 
ical modes, since reactive modes propagate always inward from the positive x side, and as 
damped outgoing modes on the negative x side. Recall that fluid modes lie asymptotically 
along direction Rj = (*,0) tr with vanishing ^-component. Taking account of this fact, to- 
gether with the faster decay rate of (0, 1)67(1, 0) tr stated in Proposition 17.11 we could by 
essentially the same argument used to prove Proposition 18 . II establish the refined result that 
the z-component of perturbation U decays faster than the u-component, reflecting the phys- 
ical picture that the fluid is in each case (weak or strong detonation or deflagration) swept 
through the traveling wave, burning completely in the high-temperature region in its wake. 
However, we do not determine precise bounds here. 
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